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Abstract. - We propose a systematic method to derive the asymptotic behaviour of the 

e^ ' persistence distribution, for a large class of stochastic processes described by a general Fokker- 

I Planck equation in one dimension. Theoretical predictions are compared to simple solvable 

I . systems and to numerical calculations. The very good agreement attests the validity of this 

y ' approach. 

o 
o 



^ Introduction. - The concept of persistence has recently motivated a lot of works, both 

p^ ' experimentally and theoretically (see Q| and references therein) : it appeared that this simple 
Cn| I concept -i.e. the probability G{t) that a random variable X{t) never goes above a certain 
CO • value (usually its mean value) during the whole time interval [0, t]- hides a real complexity, 
since its definition involves a complete knowledge of the process over a large time interval. 
Consequently, the persistence distribution gives information on the details of the process, 
^^ I which are not redundant with the traditional statistical tools, as for instance the correlation 
functions. Until now, the works was focused, besides discrete systems as spin systems 0, 
Cd I essentially on gaussian processes, and the efforts concentrated on the influence of the non 
markovian character of the process on the persistence exponent. However, it was shown 
recently p] t hat the concept of persistence can also be physically relevant for non gaussian 
situation^(^)|. Another physical situation leads naturally to study persistence in a non gaussian 
Q ] context : in nonlinear physics, the combination of discreteness and nonlinearity gives rise to 
O . localisation phenomena (the "breather- modes" ^,|[), which are often extremely pinned where 
they appear; a powerful tool to characterise the lifetime of these objects is provided by the 
persistence behaviour of the energy density on a site. In weakly coupled systems, the evolution 
of the energy density can moreover be described by a non gaussian markovian process, with 
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X: 

S . energy dependent drift and diffusion coefficients. 



(^)There is a slight ambiguity associated with the term "non gaussian" ; in P|, the process is non gaussian 
due to the non gaussianity of the noise ; but the process can be non gaussian, despite a gaussian noise, if an 
external force is present and not linear. 
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In this letter, we introduce a general method to compute the persistence behaviour for a 
large category of processes described by a one-dimensional Fokker-Planck equation. Essen- 
tially, the method deals with systems where the persistence decreases slower than an expo- 
nential law. We check the validity of our predictions by a comparison with exact results and 
numerical simulations. 

The method. - We consider X{t) a continuous markovian stochastic process. By defini- 
tion, the probability density function (pdf) ■p{x,t) = Prob{X{t) = x) obeys the Fokker-Planck 
equation 

dtP = -d,iAix)p) + ^dUBix)p) (1) 

where the functions A and B characterise the dynamical process. Once given a real value 
xo G K, the probability of (positive) persistence G+{x,t), i.e. the probability that a particle, 
which has initiated its trajectory at x > xq at time t = 0, has never crossed xo up to the time 
t, obeys the backward Fokker-Planck equation Q 

dtG+ - A{x)d,G+ + ^dl,G+ (2) 

with the boundary conditions G+{x > xq, t = 0) = I, G+{x — xq, t > 0) = 0. 

For the sake of clarity, we will first temporarily restrict the discussion to cases where 
B{x) — 2, and see later how to take into consideration the general case. In addition, we 
assume that the stationary solution of (m, Pst{x) oc expj A{x')dx' , is bounded but not 
necessarily normalisable (it means just that the potential of the force doesn't go to — oo). 

The Fokker-Planck equation (||) with constant diffusion coefficient can be mapped on a 
Schrodinger equation 0, by defining V'+(a;,t) = G+{x^t)^Jpst{x) : it reads 

9,V+ = CV'+ - (y + y)V^+. (3) 

The effective potential ^(a;) of this equation is defined on [xq, 4'0o[ by V{x > xq) — A'{x)/2 + 
A'^{x)/A and V{xo) = oo; the last precision ensures the boundary conditions ip+{xo, t) — 0, Vi 
|g|. It can be shown that the spectrum of this Schrodinger operator [—d1^ + ^{x)] is only 
located in ]0, -|-oo[ (due to the boundedness of Pst), and we know from quantum mechan- 
ics that it usually consists of a discrete part {Ai < A2 < . . . < \p} (with normahsable 
eigenf unctions V'Aj), surmounted by a continuous part ]Ac, -l-oo[ (with non normalisable eigen- 
functions V'a)- The value of Ac is usually easy to determine (since Ac is more or less related 
to liminf+oo^^(a;)/4). 

We obtain therefore the persistence distribution G+ as 

oc j 

d\ e-^*VA(x)BA \ (4) 



G+(x,t) = ^==^e-^'VA.(a;)BA,+ / 



where Bx = J^ dx ij;x{x)y'pst{x) 

The asymptotic temporal dependence of G+ can therefore exhibit three different scenar- 
ios: first, if a discrete family of eigenvalues exists for the Schrodinger operator, the tail of 
distribution is proportional to exp(— Aii), and the appropriate quantity which describes the 



Jean Farago: Persistence distributions for non Gaussian markovian processes 3 



asymptotic behaviour is the persistence time 1/Ai. Consequently, the possibihty of calcu- 
lating this time is simply related to the ability of estimating the lowest eigenvalue of the 
corresponding Schrodinger problem. 

If we consider now the cases without bound states, there are two possibilities : if Ac is 
different from zero (existence of a gap), the behaviour of G+ will not be simple a priori : 
G+{t) ex f{t) exp(— Aci) (for large t), where f{t) is unknown. An example of such a situation 
is the case of a particle experiencing a constant force toward (A < constant). In that 
simple case, we have f[t) oc t^^/"^. Our method does not deal with these cases, which we 
could call marginal because the force acting on the particle neither grows nor disappears at 
infinity. 

Finally, if A is equal to zero ("gapless" situation), we will show that there is a general 
procedure for obtaining the asymptotic behaviour of the persistence, whatever the A[x) under 
consideration. First, the formal formula 

1 f°° 

G+{x,t) = ^=== dXe-^'Mx)Bx (5) 

y^Pst[X) Jo 

shows that we have to determine the A — > behaviours of Bx and tpx in order to get the i — > oo 
limit. The cornerstone of our method is that it is possible to determine this behaviour, thanks 
to a peculiarity of the Schrodinger problem under consideration: the "A = eigenfunction" 
is known, because it is related to the stationary solution of the corresponding Fokker-Planck 
equation. The quotation are required, because this eigenfunction is highly diverging at t -^ oo 
: strictly speaking, this '0a=o does not belong to the spectrum. Before going further, let us 
show how -00 is constructed. One of the functions fulfilling ip" — {A' /2 + A^ /4:)ip = is just 
Y^pTT- Another (not proportional) is obtained by a simple quadrature, y/p7t J^ dx'p'^^{x'), and 
clearly diverges at infinity. As the "A = 0" eigenfunction must fulfil the boundary conditions 
at xq, we can write it as 

i^oix) = \/Pst{x) / — — - (6) 

JxB Pst{x') 

This "eigenfunction" is useful, since it is reasonable to assume that despite the diverging 
character of ipo, there is a continuity property of tpxix) with respect to A, which leads to the 
limit yx,\mix^Qip\ix) = V'o(2;). However, there will of course not be any related property of 
uniform convergence. 

With this assumption, which could presumably be proved using rigorous mathematical 
arguments, we have the typical portrait of ipxix) for vanishing A, over the whole range [xq, oo[ 
: "far" from the potential region, ipx{t) must become plane waves : ipxi^) ~^ Wx cos{^/Xx+4)x)- 
The phase factor (j)x is assumed to have a limit (po when A goes to zero, but its precise knowledge 
is useless for the persistence. On the contrary, it is important to know the A-dependence of 
the coefficient of proportionality Wx- The requirement that tpxix) must be normalised in the 
sense J^ dx'ipxix)iJx'{x) = (5(A - A') imposes that Wx = C* x A^^/^ §. 

On the other hand, in the region where the potential is substantially different from zero, 
the differential equation for tpx is well approximated by the equation ip" — V'tp = ; so in this 
region we can expect therefore 

^x{x)'^AxMx) (7) 

where, once again, the Ax is important to determine, by estimating the location of the 
crossover of the two limiting behaviours. This location is naturally characterised by the 
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abscissa xx given by the balance 

Vixx) ^ A (8) 

Nevertheless this condition will only hold if the computed x\ is larger than a wavelength 
~ 1/vA corresponding to the minimum typical length the eigenfunction needs to join the 
maximum of his asymptotic behaviour. In other words, if <M) gives a xx which diverges more 
slowly than 1/vA, one has to choose xx ^ 1/vA instead. In some particular cases, the 
determination of Xx is easy to achieve, since A is small. At this point xx, the continuity of ^x 
imposes that the magnitude of the two branches must be the same : this condition, 

AxM^x) ~ A-l/^ (9) 

gives the A-dependence of Ax- 

The aim is therefore achieved, because Bx can be evaluated with the approximation 

Bx^Ax dxyJpst{x)%l;o{x) + A"^/'^ / dxyjpst{x) cos(\/Aa; + (l>x). (10) 



Xq 



The relative importance of the two terms must be checked in each case. Moreover, further 
simplifications can be made ; for instance, if J y/Pst is converging, the second integral of 
( p^ ) is equivalent to 

C'' X A-1/4 / dxyf^^. (11) 



The final calculation of the persistence behaviour is achieved by the means of saddle point 
expansion of the integral 

/■OO 

G+{x,t)(x dXe-^^AxBx (12) 

Jq 

around A = 0. 

Let us consider now the general case, i.e. B{x) ^ 2. The change of variable defined by 
dy/dx = yj2/ B{x) = g leads to a new Fokker-Planck equation for the distribution p{y) = 
p{x)dy/dx, with the new coefficients B{y)/2 = 1 and A{y) — A{y)g{y) + dyg/g. As a function 
of y, we have then to compute the characteristic abscissa yx- It is however more convenient 
to compute the corresponding abscissa in the original coordinates, i.e. the xx defined as the 
solution of 



1 B(x\) dA , , 1 c-Q . s , , 



where A{x) = J^ x (^A{x) ~ \d^B{x)\^ . (14) 

The "00 eigenfunction can also be expressed in the original x coordinate as 

/ : r^ OAq.' 
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With this change of variables, the coefhcient Ax is given by Ax^o{xx) ^ A^^^'*. The same 
translation can be performed on Bx, leading to 



Bx=Ax dx{2/B{x)Y'^^p,t{x)Mx) 



+ C' X A-i/4 / dx{2lB{x)f'^^p,t{x) cos{VXy{x) + 0o) (16) 



The presence of y{x) in the cosine is not a problem in general, because only its asymptotic 
behaviour is needed, which is often easily computed. 

In summary, the cases where the diffusion coefficient B{x) depends on x are not a handicap 
to the procedure proposed in this letter. 

Examples. ~ To confirm the above procedure, we will first consider exactly solvable 
examples. The simplest is probably the free brownian motion, i.e. A = 0,B/2 = 1. In that 
case, V = and xx behaves as 1/VA (the prescription (g) giving xx = 0, we have to choose the 
"minimal divergence" I/VA) ; moreover, ^o{x) "^ x, a.nd pst = C**. Consequently, Ax ^ X^^'^, 
Bx ~ A"^''"'^ and the persistence goes like t~^/^, as the exact result states. 

To check the theory on a less trivial situation, let us consider the cases B/2 — 1,A = 
—vx^^ with V > Q. The Schrodinger potential is proportional to a;~^, and it is easy to 
verify that, for large x, Pst oc a;"'', -00 « W^+^xa cx \-^I^,Ax cx X^^^+^^'^^Bx cx A^'^-^)/''. 
This gives the result G+{x,t) cx t~(i+'^)/2. The Schrodinger problem can actually be solved 
exactly, because the eigenfunctions ^px are proportional to \fx{J ^yj^iy2{\Xx)Y(^yj^iy2{\Xxa) — 
J{v+i)/2{vXxQ)Y/^j^iy2{'yXx)) llfl]. The exact result matches our prediction (and the limit 
V ^r Q gives the free diffusion result). 

We turn now to situations for which the knowledge of the spectrum is unknown. For 
instance, consider the cases A — —vx^"' , B/2 — 1, with i^ > 0, a > 0. 

If we consider first the cases a > I, we have then (for large x) Pst oc exp(i/x^^"/(a — 1)), 
V'o (X a:, y cx a;-"-\ xx cx X'^/^°'+^^ ^ xa cx X'^/"^, Ax cx X^l^, Bx cx A-3/4. it gives 
G{x,t) ^ i^^/^. It is interesting to note that the potential slope is inefficient to modify the 
persistence behaviour of the particle, and that one obtains a discontinuity of the exponent as 
a function of a (when a —^ 1). 

For cases < a < 1, the situation is completely different : ipo{x) ^ a;" exp(i/a;^~"/2(I — 
a)), V - a;-2" and xx - X'^^^°'. It leads to Ax ~ Ai/4exp(-(A/Ao)-(i-")/2a)^ ^^ _ 

AaA-(i+""')/^ where Aq cx 4(i//2(f - a))^'^-"'>/^" /u^. Finally, the persistence is found to 
behave like 

G ^ exph(i/io)('-"^/('+")] X t-(3"-i)/2(a+i) 

with to cx ,y2/(i+a)4(a-i)/(a+i)/(;^ - q) X [C^" + C""^] (( = (1 - a)/2a). The coefficient of 
proportionality for tg must be of order I. 

To check the validity of our theory, we have numerically computed the persistence distri- 
bution for the last case (a < 1) : figure |l| (a) shows three different curves of log(G') versus 
Vi for a = 1/3 (note that in that case, the persistence is purely a stretched exponential, 
with exponent 1/2) and different values of v ; the inset shows the v dependence of the slopes 
of these curves, which is well predicted by the expression of tg given by the above formula. 
Panel (b) of the figure shows two others examples of log(i^^"~^-'/^^""''^'G') plotted as function 
of C"** X t(i~")/(i+")j emphasizing again that the linear behaviour analytically derived is fully 
satisfied. 
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Fig. 1 - Persistence distribution numerically computed for A{x) — —ux'" with < a < 1. (a) 
a = 1/3 and v — 0.5, 1, 2. The inset shows the ^/-dependence of the coefficient to (stars) compared 
to the prediction 1/'^''-^+°'' of the theory (circles), (b) a = 0.5 (dots) and a = 0.75 (circles) (in that 
panel, the values of to are of order 0.12). The ultimate diverging tails of these curves are artefacts 
due to a poor statistics in these regions. 



As a last example, let us consider now the case of an underdamped particle in a potential 
well, submitted to a damping 'y{E), a priori function of its energy. Kramers |ll|,|l^ had shown 
that the fast angle variable can be eliminated, leading to an effective Fokker-Planck equation 
for its energy E, with A = ^[kBTdEh{E)I{E))-j{E)I{E)] and B = 2kBT^-f{E)I{E) 
(where I{E), uj{E) are the action and the pulsation of the trajectory). If the damping vanishes 
rapidly enough as E increases, the problem belongs to the "gapless" category. For instance, 
if one considers the case of an harmonic potential {ui = ujq, I{E) = 2ttE/luo), with a damping 
'y{E) = E~°' (for E sufficiently high), one has that, if a > 1, the Schrodinger potential goes to 
zero at infinity, and the persistence is asymptotically G ^ /;3/2q-i exp — (i/ig)" ■ This class of 
situations is physically particularly relevant to the cases of non linear coupled oscillators: the 
frequency shift between two adjacent oscillators having different energies leads to a enormous 
slowing down of the diffusion of the energy jlj, Q , described in a mean-field model by such 
a vanishing damping. 

Correlation functions. - It is interesting to remark that this method can be applied to the 
calculation of tails of correlation functions, for systems belonging to the appropriate "gapless" 
class. The hypothesis that a correlation function exists implies that a real equilibrium is 
reachable by the system ; with this restriction, the correlation function of the variables is [|6| 



< x{t)x{0) > = / dX 



'\t 



dx x\/pstix)ip\{x) 



(17) 



It is clear that an analogous derivation can be performed on that formula, in order to extract 
its asymptotic behaviour with, nevertheless, a slight difference in the definition of the i/'a : 
the boundary condition ip\{xo) — does no longer exist, and the "00 function is now a real 
eigenfunction equal to y/pit- 
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Conclusion. - In this letter, we have presented a general method to derive analytically 
the asymptotic behaviour of the persistence probability, for a large number of markovian 
processes in one dimension : we have shown that it is possible to classify markovian processes 
in three categories ; the first is characterised by an exponential extinction of the persistence, 
a second one, quite marginal, whose treatment is beyond the scope of this letter, and the 
third, which we called "gapless" , in reference to the structure of the spectrum of an associated 
Schrodinger operator. For this third class, we have introduced a systematic procedure to 
obtain the persistence distribution at large times, and we have tested the validity of the 
procedure by two complementary ways : we compared the results of our method to exactly 
solvable models, and the theoretical predictions to numerical simulations, when the exact 
result is not known. Both comparisons have shown an impressive agreement. 

It is interesting to note that this method could be presumably extended to the treatment of 
multidimensional cases, as soon as the corresponding Fokker-Planck equation can be mapped 
with an appropriate change of variables on a Schrodinger equation. In this extended version 
of the theory, relevant quantities like xx would become functions of the solid angle of the 
parameter space, leading probably to a more complicated behaviour of the persistence. This 
extension could therefore be an original way to study non markovian cases, since memory 
effects can always be interpreted as an elimination of additional variables. 



I thank T.Dauxois and A.Alastuey for comments and discussions. 
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